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Abstract 

We have theoretically investigated the spin and charge fluctuations in the 
quasi-one dimensional organic superconductor (TMTSF)2C104. Using the 
extended multi-site Hubbard model, which contains four sites in a unit cell 
and the transfer energies obtained by the extended Hiickel method, we cal- 
culate the linearized gap equation with the random phase approximation, 
to flnd novel order parameters of superconductivity due to several kinds of 
charge fluctuations induced by the anisotropic intersite repulsive interactions. 
For the singlet state, the order parameter with line nodes appears in the case 
of the strong charge fluctuation, while the order parameter with anisotropic 
gap suggested by Shimahara is reproduced in the spin fluctuation. The triplet 
state is also obtained for the wide parameter range of repulsive interactions 
due to a cooperation between charge and spin fluctuations. 

Keywords: Organic Conductor, TMTSF, Superconductivity, Electron 
Correlation, Spin Fluctuation, Charge Fluctuation 



1. Introduction 

In quasi-one dimensional organic conductors TMTSF (tetrametyltetrase- 
lenafulvalence) salts, there are several experimental evidences of supercon- 
ductivity, [H, i, 3 i, i, B 0, B B [io|L which exhibit the large critical held being 
much larger than the Pauli limitll, B B [13 ^or (TMTSF) gPFg, the triplet 
superconducting(SC) state has been maintained from the measurement of 
the Knight shift and the NMR relaxation rate. 0, For (TMTSF)2C104, 
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[l[ in the weak magnetic field, the singlet state with line nodes has been sug- 
gested by the the Knight shift and relaxation rate of the NMR experiment in 
the SC phase. 0,0] Further, it has been suggested that the triplet or FFLO 
state emerges under strong magnetic field, jjlsl, 0, |lo[ 

The specific property of (TMTSF)2C104 emerges due to the anion or- 
dering below 24K. The folded Fermi surfaces with a unit cell containing 
four TMTSF molecules is obtained from the transfer integrals, which are 
estimated from the extended Hiickel method based on the X-ray diffraction 



measurement. 



11, 12, 13 



Several theoretical works have been performed to understand the super- 
conductivity in the quasi-one dimensional organic conductors. The single 
band models for (TMTSF^PFg has been investigated by perturbation the- 



ories. 



]41l5|,ll6 



J ., 17, isl. 19| The magnetic field favors the triplet or FFLO 

superconductivity. 

Ha 

llTI , |18| while the on-site repulsion induces the singlet 
superconductivity in the absence of magnetic field. [l9[ The interplay be- 



tween the spin and charge fluctuations in the triplet superconductivity has 
been investigated by the renormalization group theories. 20|, |2l|, |22| The 



inter-site repulsions play important roles for the triplet superconductivity. 
3 E [21I Q In the muhi band models for (TMTSF)2C104|23l, [3, on the 



other hand, the singlet superconductivity with anisotropic gap (but without 
line nodes) has been suggested in the absence of magnetic field. [2^ 

In the present study, the possible singlet state with line nodes in the 
absence of magnetic field is investigated, since the experiment indicates the 
existence of line nodes, jsj We use the extended multi-site Hubbard model rep- 
resenting the system of (TMTSF)2C104, when the transfer energies are given 



by the extended Hiickel method with the X-ray diffraction measurement. [13 



The unit cell consists of four TMTSF sites owing to the presence of the 
anion ordering. Fruther, the roles of the charge fluctuation for the triplet 
state are also investigated on such a multi-site system. The spin and charge 
fluctuations are treated by the random phase approximation (RPA) on the 
site-representation, and the order parameters of the superconductivity 
are evaluated using the linearized gap equation as described in §2. The singlet 
and triplet superconductivities are investigated by varying the magnitude of 
the anisotropic nearest-neighbor repulsive interactions with the fixed on-site 
interaction in §3. The origins and properties of those SC states are discussed 
in §4. In §5, summary and discussion are given. 
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Figure 1: (a) A model describing the two-dimensional electronic system for 
(TMTSF)2C104. The x-axis is taken as the stacking direction. The unit cell (the dotted 
line) consists of 4 molecules with several transfer energies id;a/3- (b) Inter-site repulsive 
interactions Va,Vb and V , where V is discarded. 



2. Formulation 

2.1. Model 

Figure [1] (a) sliows tlie transfer integrals describing the two-dimensional 
(2D) electronic system for (TMTSF)2C104, where we choose the x-axis along 
the stacking direction of TMTSF molecules. The thickness of the lines con- 
necting between the molecules represents the relative magnitude of td;a/3- 
Figure [1] (b) displays the inter-site replusive interactions, which are taken in 
the present model. 

Based on Figs. [H (a) and (b), the extended multi-site Hubbard model is 
given by 

n = no + n', (i) 

with 

^0 = ((^d;a/3 - /i'5a,/3'5(d))aL<7«i+d/3<7 + h.c), (2) 

ido/3(T 

and 



"H' - ^ f/a^a^ai„4.aiat + ^ V'd;a/3aL(7«i+d/3(7''2i+d/3<7'aia,7, (3) 

ia idal3acr' 
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where i(= (1, 1) ■ ■ ■ {y/Ni, y/N^)) denotes a two dimensional lattice vector 
for the unit cell, and d is a vector connecting different cells. Quantities 
a and /3(= 1,2,3,4) are indices of molecules in the unit cell and a{=^,l) 
represent spin. The transfer energies, td-a/s, are obtained by the extended 



Hiickel method based on the X-ray diffraction measurement. 1 121. 1 131. 1 111] Cou- 
pling constants U and Vd-ai3 correspond to the on-site and nearest neighbor 
Coulomb interactions, respectively. For Vd-a^, there is a relation, Vd-^a/s = 
= V-d;ai3- The uuit of the energy is taken as eV. 
We investigate the superconductivity by varying the nearest neighbor 
Coulomb interactions between the TMTSF chains, i.e., Va along the y axis 
and Vb along the x + y axis (see Fig. [I](b)). In the RPA, the critical value of 
the nearest neighbor Coulomb interaction along the intrachain, V, at which 
the charge fluctuation diverges, is more than twice of those of Va and H- 
This comes from the fact that the role of the nesting vector is incompatible 
with that of V'. The same fact has been shown in the renormalization group 



theory. [22| Thus V is discarded in the present study. 



Using Fourier transformation, the Hamiltonian is rewritten as 

kcra/3 

+ ]^ f^ajc+qta^k'-qia^k'^aakta 
kk'qcj 

+ OAT X] ^a/3(q)'^k+qaQ'^k'-qa'/3«k'<7'/3akaa, (4) 



2N 

kk'qccr'Q, 



where 



ei2(k) = ti2 + t[2e 
ei4(k) = hA + t\^e-'^\ 



£23 (k) = t23 + 4e"'''^ 

e24(k) = t24 + 4(e"*'^ + e-''^)+4e"*^'^^'^\ 

e34(k) = t34 + t'34e-*'^ (5) 



K3(q) = Vu + Vl^ie'"^- + e-"'y) + V;'^e'^''--''y\ 
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^23(q) 
V^24(q) 



+ V'e 



23 -r »'23 
^24 + Vi^{ 



■iqy 



-tqx 



+ e 



-iqy- 



+ v'' 



24^ 



-i(qx+qy) 



with e/3o 

energies are given by 1 121. 1 131 . Ill 



V*j^ and the other elements are zero. 



t 



12 



t' 



34 



0.413, t'l 

0.362 along conduction chains, ^23 = ^' 



"14 



12 - 0.324, t34 

= -0.050, ti4 = 



along y-axis, and t'24 = 0.070,^24 = ^" 



24 



0.020, t 



13 



13 



(6) 

The transfer 
= 0.335 and 
t'23 = -0.100 

= 0.071, t'lg = 



0.021 where t, t' and t" correspond to transfer energies between molecules in 
a unit cell, those between the nearest neighbor cells, and those between the 
next-nearest neighbor cells, respectively. Note tai^ = tpa- We choose V14 = 
V23 = V^i4 = ^3 = Va for nearest neighbor, i,e, d = (±1, 0), (0, ±1) and V13 = 
V24 = V24 = V/s = Vft for next nearest neighbor, i.e., d = (±1, ±1), (±1, =f1). 
(See Fig. [D^b).) The tight-binding Hamiltonian Eq. ^ is diagonalized by 



J](e„^(k)-/iM<(k) = eX(k), 



(7) 



where A{= 1, 2, 3, 4) is an index of the band representation and ^"^(k) denotes 
the eigen energy measured from the Fermi energy (= /i). note that the 
eigenvector c?^(k) gives the transformation from the site representation to 
the band representation. 

In Fig. El^a), the energy band spectrum is shown. There are four bands 
due to four molecules in the unit cell where the first and second bands from 
the top crosses the Fermi surface. The shape of quasi one- dimensional Fermi 
surface is shown in Fig. [2Ib) representing the folded Fermi surface. In Fig. 
131 the density of states ,-D(e), is shown where 



Die) 



Nti ^ \t 



1 



16 -^^ + 11 



(8) 



and sum rule is given by J D{t)dt = 8. These double one- dimensional Fermi 
surfaces have been observed by AMRO[ll|. 



2.2. Pairing Interactions and Gap Equation 

The spin and char ge fl uctuations in the weak coupling regime are treated 



by RPA. |14| . Il5l . Il6l . Il7l Il9l . l26l |25[ where the site representation [26j for 



singlet and triplet pairing interactions are given respectively as 

pSinglet ^fj^y^ ^Ux^U - ^{U + 2V)x^\U + 2V), 



(9) 
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pTViplet _ \lJj^S(j _ hfj + 2V)x^{U + 2V). (10) 

The hat denotes a matrix representation with an element (/)«/? = /o^g. The 
pairing interactions are obtained from -y^ and x*" which denote spin suscep- 
tibihty and charge susceptibihty; 

x'-{i-ftj)-^x\ (11) 

and 

X^ = (7 + 2x°T> + x°J7)-Y. (12) 
The quantity x° is an irreducible susceptibility which is given by 

(x°(q))a/3 - ^5^d^(k + q)d;^(k + q)C(k)df(k) 

/(Ck+q) ~ /(Ck) /,o\ 
Sk+q Sk 

where /(x) = (e^/^ + 1)~^ is the Fermi distribution and T denotes tempera- 
ture. 

In order to examine a symmetry of order parameter of the superconduc- 
tivity, we solve numerically a linearized gap equation given by, 



ASa/,(k) = --^ [(^^(k - k')).;.GL,(k')G?^,(-k')] S^^(k'). (14) 



k'a'/3' 



E^^(k) is anomalous self energy and G°Q,,(k') denotes a bare Green's function 
written as 

G'°^(k,i£n) = E^a(kX(k)-^, (15) 

where = (2n + l)7rT is Matsubara-frequency {ks = 1) and n is an integer. 
The quantity (A = Singlet, Triplet) denotes the paring interactions for 
singlet state and triplet state which are rewritten as 

pSinglet _ pSs _|_ pC ^-|^g^ 

and 

pTriplet _ pSr _|_ pC ^-^j^ 
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The r.h.s. of Eqs. (fTB]) and (fT7|) is given by 

P^' = tj+^Ux'U, (18) 

P'^ = -lux'^U, (19) 
P^ = V -1-{U + 2V)X^{U + 2V). (20) 

Note that P"^^ and P*^^ represent spin fluctuation, and P*" denotes the 
charge fluctuation. 

The transition temperature for the SC state, Tc, is determined by the 
condition, A = 1. The gap function at T = 0, Ao^(k), is estimated from 
Ao^(k) = 5fE^^(k), where g is chosen to reproduce the gap which is ex- 
pected from Tc. ji^ Using AQ,^(k), we get the quasi-particle bands by 
diagonahzing Hamiltonian 



where 



^ - eQ:/3(k)akQ,^ak/3cr 

ka/3(T 

- ^ A„^(k) a^a^_k^^ 

~ X]^^/3*^'^)"-k/3iakat' (21) 

ka/3 

X^e„,,,(k)<(k) = Eif'd5(k), (22) 

/3' 



_ ( e(k) - - A(k) , . 

'^''^"^ -At(-k) -*e(-k)+M.^ ^ ^ 



3. Singlet SC state vs. triplet SC state 

3.1. Band Calculation and Model Setting 

In the present study temperature is fixed at T = 0.01 [eV]. We take 24 x 
24 meshes for calculating the bands, the susceptibilities, and the supercon- 
ductivity. The energy bands are shown in Fig. [2](a). Since the electrons are 
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Figure 2: (a) Energy Band for (TMTSF)2C104 given by Eq. © with the 3/4 filling, 
where the two bands 1 and 2 from the top cross the Fermi energy, /i. (b) The Fermi 
surfaces for the first and second bands, where the arrow represents the nesting vector. 



3/4- filled, the Fermi surfaces exist in upper two bands. The density of states 
(DOS) is shown in Fig. [3l There is a large Van Hove singularity above the 
Fermi energy. 

Now we consider the electron correlation effects. It is expected that 
U plays important roles for superconductivity in (TMTSF)2C104, because 
(TMTSF)2C104 is located next to the spin density wave (SDW) state in the 
phase diagram suggested by Jerome. Then we calculate the possible SC 
states in (TMTSF)2C104 by choosing U = 0.60 which is close to the critical 
value for the divergence of the spin fluctuation in the RPA. We investigate 
the roles of the charge fluctuation by varying the interchain repulsive inter- 
actions, Va and Vb, with the fixed U, where we define V as 

V ^ + (24) 



3.2. Superconductivity for V =0 

For Va = Vb = 0, spin and charge pairing interactions {P^ and P'") are 
shown in Fig. H] where is much larger than P'-^ . The susceptibility 
which determines P^ is strongly enhanced by U compared with x*^, since the 
q vector for the peak of is the same as the nesting vector. 

Now we examine numerically Eq. ( !T^ to find that A for spin singlet is 
larger than that for spin triplet, e.g., A = 0.9 for the singlet state and 0.3 
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Figure 3: The density of states (horizontal axis). 



for the triplet state at U = 0.6 and K = H = 0. In the present paper, a 
phase factor in anomalous self energy is determined in order to maximize the 
value ^a/3k l-^^^a/?^^)!- "^^^ anomalous self energy consists of 16 elements 
of {a, P) on the basis of site representation. Diagonal elements {a, a) are 
larger than others. The momentum dependence of these elements is given 
by "coskx" which indicates a sign change on a line separating two Fermi 
surfaces in Fig. MJo). The absence of the node on the Fermi surface, which 
gives a full gap state ( written as "singlet (full)"), is consistent with that of 
Shimahara's result. |23'] As for the off-diagonal element of E^^(k), there is a 
clear contribution from the same chain , e.g., (1,2) and (3,4), while that from 
the interchain is negligibly small. These behavior of intrachain dependence 
suggests a role of dimerization, which reduces to the half-filled band as shown 
in Fig. H 

3.3. Superconductivity for V ^ 

We examine the effect of inter chain interactions V to find a possibility 
of new SC state which is induced by charge fluctuation. With increasing V, 
X*^ increases, and diverges ai V = Vc- In order to calculate the region where 
the charge fluctuation becomes a dominant one, V is chosen to satisfy the 
condition that x*" > Sx"^ and V < Vc- We examine the SC state on the line 



from {Va,Vb) = {0,0) (O) to several points given by (K, Vb) = (0.00,0.64) (a), 
(0.11,0.64) (b), (0.21,0.64) (c), (0.33,0.62) (d), (0.49,0.49) (e), (0.67,0.22) (f) 
and (0.70,0.00) (g), respectively. The ratio of \4 to VJ, increases as the point 
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Figure 4: The momentum dependence of pairing interactions for the elements Pj^ and 
P^i at T = 0.01 and ^ = 0. The red and blue represent the repulsive and attractive 
interactions, respectively. The symbol O corresponds to Va = Vb = 0. 



moves from (a) to (g). By solving Eq. (fT4l) . A for singlet and triplet states 
are calculated and the larger one is chosen as the expected state. 

In the region of small V, the singlet SC state similar to one at = is 
reahzed where the state with a full gap is found. In Fig. [5l the V dependence 
of A of both singlet and triplet states is shown as a function of V for three 
cases of (O) — )■ (a), (O) — )■ (d), and (O) — )■ (g). For the variation from (O) 
to (a), (the first case) , the singlet state is suppressed and the triplet state 
is expected. For the variation of V from (O) to (d) (the second case), the 
singlet state moves to the triplet state close to = Vc-. For the variation of 
V from (O) to (g) (the third case), the singlet state remains as the dominant 
state but moves to another singlet state (as explained later) after A takes a 
minimum close to (g). 

Here we note the V dependence of A for triplet state, which is enhanced 
by x'-' . From Eq. ffTOj) . the term for the spin fluctuation and that for the 
charge fluctuation has the same sign in triplet pairing interaction. Although 
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Figure 5: V dependence of A for singlet and triplet states for T = 0.01, C/ = 0.6. The 
horizontal axis is shown for V/Vc > 0.7 which represents from the point (O) to (a), from 
the point (O) to (d) and from the point (O) to (g), where these points in the parameter 
space are given by (Va, Vf,) = (0,0) (O), (0.00,0.64) (a), (0.33,0.62) (d), and (0.70,0.00) (g), 
respectively. The solid line denotes A for the singlet state and the dashed line denotes A 
for the triplet state. 



the triplet SC state is enhanced but the singlet state is suppressed by charge 
fluctuations, such an effect becomes large only close to ^ = Vc- This is why 
the singlet state is favored and the region for the triplet state is small even 
in the presence of V. 

4. Pairing interaction and Anomalous self-energy 

Now we examine the respective SC state precisely based on the pairing 
interaction shown in Figs. |6] and anomalous self-energy. 

(i) First, we discuss the case indicating a possible triplet SC state, which 
is given by (a) and (d). When x'^' is increased by V, the singlet state is 
suppressed but the triplet state is enhanced (Fig. [6] (a) and (d)). Actually, 
the competition between these two kinds of fluctuations occurs for the singlet 
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Figure 6: The momentum dependence of charge-fluctuation induced pairing interaction 
(Pii) for V ^ 0. where parameters corresponding to points (a), (d), (e) and (g) are 
(0.00,0.64), (0.33,0.62), (0.49,0.49) , and (0.70,0.00) , respectively. The attractive interac- 
tion is represented by blue on the kx-ky plane. For the comparison, P^^ is shown. 



state as seen from Eq. (Q and by noting a fact that the peak position of 
P*-^ is nearly the same as that of Pf from Fig. HI The triplet state may 
be expected when V ^ Vc, and is rather favored by the introduction of Va. 
When the point moves from (a) to (d) (i.e., increasing Va), the peak position 
of P*" begins to deviate slightly from (vr, vr) and the attractive interaction 
increases due to the increase of the amount of P*" integrated with respect to 
the momentum. 

(ii) Secondly, we study the region where the singlet state is favored by 
the effects of x*" (points (e) - (g). The pairing interactions (P*") induced by 
charge fluctuation are shown in Figs. |6] (e) and (g) where the position of the 
peak is different from that of in . As a result, the effect of x'" on the singlet 
pairing state becomes compatible with that of x^ for V ^ Vc- Further since 
P*-^ is larger than P^, the spin singlet state is determined mainly by x'" and 
is assisted by x^ . The change of the self-energy along axis is very small. 
This comes from the difference in the peak position of charge fluctuations 
i.e., {n,n) for (a) and that on the axis for (e) and (g). 

Finally, we examine the node of the gap, which is calculated from Eqs. 
f l2T|) and f l22|) . Figure [7] denotes the position of nodes and Fermi surfaces for 
points (O), (f), (g) and (d). For the singlet state (O), which is mediated by 
X^ the node exists between two different Fermi surfaces implying a full gap 
due to the existence of the two bands. j23| For (f) and (g), the singlet state 
mediated by x*" has two nodes at ky = ±vr/2, which cross the Fermi surfaces. 
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Figure 7: Fermi surface and gap node where the dashed hues denotes Fermi surface and 
the sohd hne denote gap node at T = 0.01, U = 0.6. 



Gap function for triplet has a node on the hne = suggesting a fuU-gap 
state. 



5. Summary and Discussion 



We have found the possible singlet SC state with the line nodes in the 
extended multi-site Hubbard model representing (TMTSF)2C104 with the 
multi-band and under the anion ordering. Such a state is induced by the 
charge fluctuation owing to the interchain Coulomb interaction Va (along the 
y-axis), while the singlet SC state with anisotropic (but finite) gap suggested 
by Shimahara 23|] is reproduced in the presence of the strong spin fluctuation. 
Moreover, we have investigated the roles of the charge fluctuation for the 
triplet SC state. It is found that the triplet state appears in a certain range 
of the interchain Coulomb interactions, Va and Vb, It is noted that the charge 
fluctuation cooperates with the spin fluctuation for the triplet pairing, but 
competes for the singlet pairing. For such a triplet state, the p-wave like state 
is obtained in the present multi-site model, while the f-wave gap function have 



IJ, ll5|, ll6|, ll7|, ll8|, |20|, l2l|, |22 



been suggested in the single site models. 

Finally we discuss how to distinguish the SC state with the lines and 
that with full gap of anisotropy on the Fermi surface. Using a maximum, 
Amax5 and a minimum, A^i^, for the latter case, 1/Ti reveals the exponen- 
tial behavior at low temperature, i.e., T < Amin- In the preset case, it is 
expected that 1/Ti exhibits exponential decrease at lower temperatures than 
TcAjnin/Amax ~ 0.2K uslug the ratio given by both our and Shimahara's 
results. On the other hand, Yonezawa et al. have reported that the 1/Ti 
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exhibits behavior at temperatures lower than 0.2K[8|. This fact suggests 
that the SC state with the hue nodes supports the result of the NMR exper- 
iment in the absence of magnetic field. 
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